arXiv:astro-ph/0112254vl 11 Dec 2001 


High-Resolution Spectroscopy of G191-B2B in the Extreme 

Ultraviolet 


R. G. Cruddace\ M . P. Kowalski^, D. Yentis^ C. M. Brown^ H. Gursky^ 

M. A. Barstow^, N. P. Bannister^, G. W. Fraser^, J. E. Spragg^ 

J. S. Lapington^, J. A. Tandy^, B. Sanderson^, J. L. Gulhane^ 

T. W. Barbee^ J. F. Kordas^ W. Goldstein^ 
and G. G. Fritz® 

ABSTRACT 

We report a high-resolution (R=3000-4000) spectroscopic observation of the 
DA white dwarf G191-B2B in the extreme ultraviolet band 220 - 245 A. A low- 
density ionised He component is clearly present along the line-of-sight, which if 
completely interstellar implies a He ionisation fraction considerably higher than 
is typical of the local interstellar medium. However, some of this material may 
be associated with circumstellar gas, which has been detected by analysis of the 
GIV absorption line doublet in an HST STIS spectrum. A stellar atmosphere 
model assuming a uniform element distribution yields a best £t to the data which 
includes a significant abundance of photospheric He. The 99-percent conhdence 
contour for the fit parameters excludes solutions in which photospheric He is 
absent, but this result needs to be tested using models allowing abundance gra¬ 
dients. 
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1. INTRODUCTION 


White dwarfs are among the oldest objects in the Galaxy. As remnants of all stars 
with an initial mass < SM©, they are important laboratories for study of evolutionary 
processes and the behavior of matter at extreme temperature and density. Study of their 
space and luminosity distributions helps map the history of star formation and could, in 
principle, determine the age of the disk, yielding an important lower limit to the age of the 
Universe. Further, cool white dwarfs may account for a substantial fraction of the missing 
mass in the galactic halo (Oppenheimer et al. 2001). However, these goals depend on our 
understanding of white dwarf evolution, and in particular on predictions of the cooling rates. 
These in turn are affected by the mass, radius and photospheric composition of the stars. 
Our understanding of the physical mechanisms that determine white dwarf evolution leaves 
unanswered several major questions. While the emergence from the AGB of two groups of 
white dwarfs, whose compositions are dominated by H or He, is beginning to be understood, 
the complex relationship between these branches, and a demonstrable temperature gap in 
the cooling sequence of the He-rich branch, cannot yet be explained. Determination of the 
photospheric He and heavy element content would provide important information on the 
evolutionary history of the stars. Already it has been established that significant quantities 
of elements heavier than He are present in the atmospheres of the hottest (T > 50, 000 K) 
white dwarfs (Barstow et al. 1998). 

G191-B2B is one of the brightest and best-studied of the hot H-rich DA white dwarfs. As 
it lies near the top of the DA cooling sequence, measurements of its effective temperature, 
surface gravity and composition represent an important benchmark in the study of the 
whole DA sample. lUE and HST observations at high spectral resolution have made it 
clear, that G191-B2B falls into the group of very hot DA stars (T > 50,000 K) whose 
atmospheres contain signihcant quantities of heavy elements. In particular, detections of 
G, N, O, Si, P, S, Fe and Ni have been reported in various studies (Bruhweiler & Hondo 
1983; Sion et al. 1992; Vennes et al. 1992, 1996; Holberg et al. 1994). Such material is 
responsible for severe depression of the extreme ultraviolet (EUV) flux in G191-B2B at 
A < 200 A, when it is compared with stars having pure H atmospheres, and the star has 
been an important target for spectroscopic EUV observations to determine the principal 
opacity sources. In addition, an important goal is to obtain a self-consistent model, having 
an effective temperature, surface gravity and composition which can £t the far ultraviolet 
(FUV) and EUV observations simultaneously. This would demonstrate our understanding 
of the star and the reliability of the model calculations, which could then be applied to other 
objects. 

However, a complete understanding of the EUVE spectrum of G191-B2B has been 



3 


elusive. Initial attempts to match the observation with synthetic spectra failed to reproduce 
either the flux level or the general shape of the continuum (see Barstow et ah (1996)), as the 
model contained insufficient Fe and Ni lines. Consequently ~ 9 million predicted lines were 
added to the few thousand with measured wavelengths, yielding a self-consistent model able 
to reproduce the EUV, FUV and optical spectra (Lanz et al. 1996). However, good agreement 
could be achieved only by inclnding a signihcant quantity of He, either in the photosphere or 
in an ionised interstellar component. Unfortnnately, due to the limited resolution of EUVE 
(~ 0.5 A) in the He H Lyman series band, the He contribution could not be detected directly. 
More recently, it has been shown that photospheric heavy elements may not be distribnted 
homogeneously in the radial direction, so that more complex stratihed structnres should be 
considered. This approach has yielded a good £t of model atmospheres to observational data 
across the soft X-ray, EUV and FUV bands (Barstow, Hnbeny & Holberg 1999; Dreizler & 
Wolff 1999). Important progress has been made also in incorporating radiative levitation 
and diffusion self-consistently into the calculations (Dreizler & Wolff 1999; Schuh, Dreizler 
& Wolff 2001). The need for a He contribution is rednced in these stratified models but 
not eliminated. We present the results of a search for the He H component in the EUV 
spectrum of G191-B2B, using the J-PEX (Joint Astrophysical Plasmadynamic Experiment) 
high-resolntion EUV spectrometer. This spectrnm was obtained by J-PEX after lannch by 
a sounding rocket in Febrnary 2001. 


2. The J-PEX High Resolution EUV Spectrometer 

The J-PEX spectrometer, described by Cruddace et al. (in preparation), is a slitless, 
normal-incidence instrument employing hgnred spherical gratings in a Wadsworth mount. 
The optic comprises four ion-etched laminar gratings with a groove density of 3600 gmm”^ 
and a focal length of 2.0 m. The gratings are coated with a MosC/Si/MoSB multilayer 
designed to operate in the band 220 — 245 A, and are optimised to suppress zero order and 
to yield maximum efficiency in first order at 235 A. At this wavelength the spectrometer 
achieves an effective area of 3.0 cm^. The design specification for the spectral resolving 
power (R) in flight is 4980, which includes the effect of a pointing uncertainty of 1 arcsec. 
Calibration of R was hindered by thermal deflections in the spectrometer, caused by detector 
heating over long exposnre times, and consequenctly the average measured resolving power 
was 2750, yielding an estimate for inflight resolving power of 2600. However the problem 
encountered in calibration should have a much smaller effect during a flight lasting only 
5 minutes. Therefore we can only set reasonable bonnds to the resolving power in flight, 
namely 3000 < R < 4000. This uncertainty is taken into account in analysis of the flight 
data. J-PEX was launched by a NASA Black Brant IX sonnding rocket (NASA 36.195DG) 



at 05.45UT on 22 February 2001. The payload completed its mission successfully, during 
which it observed the target G191-B2B for 300 seconds. 


3. Data Reduction and Analysis 

The spectra of the four gratings were recorded independently by the focal-plane detec¬ 
tor. In Cruddace et al. (in preparation) we describe how the positions of photon events in 
each spectrum were corrected for ACS drift and jitter, and the wavelength of each spectrum 
was calibrated. The events were summed in bins of width 0.024 A, before being superposed 
to yield one spectrum. The bin width was chosen to oversample the data by a factor 2.4 in 
comparison with the resolution at R = 4000, so as to minimise loss of information during 
superposition. The hnal spectrum, in which the signal-to-noise ratio (S/N) has been max¬ 
imised by increasing the bin size to 0.048 A, is shown in Figure 1. The average S/N per 
bin was 5.0. The background count-rate in the detector field during the observation was 
4.2cts“^, sufficiently low that the spectra were essentially free of background. Therefore 
the error bars in Figure 1 have been assigned assuming Poisson statistics. The spectrum 
contains one flaw caused by the pointing on target, in which data above ~ 239 A was lost 
in two of the spectra. Thus only half the instrument effective area was used in this region, 
in which accordingly we have increased the bin width to 0.096 A. 

We have compared the observed spectrum with the predictions of a model, in which the 
white dwarf atmosphere has a homogeneous composition and the flux is absorbed by H I, 
He I and He H in the interstellar medium (ISM). Although stratihed atmosphere models are 
more successful in reconciling spectra of G191-B2B in the EUV and FUV bands (Barstow 
and Hubeny 1998; Barstow, Hubeny & Holberg 1999), the homogeneous models are adequate 
for the relatively narrow J-PEX waveband, and give a useful baseline for comparison with 
abundance measurements made in other studies. The analysis technique has been described 
extensively in earlier papers (e.g. Lanz et al. 1996), and we give here only a brief overview. 
The XSPEG software was used to fold model spectra through the J-PEX instrument response, 
and as we were dealing with a spectrum having a small number of counts per bin, the 
best match between model and data was obtained by minimisation of the Gash statistic 
(Gash 1979). This does not assign an absolute value to the goodness of fit, but does allow 
uncertainty ranges to be determined for each free parameter. 

The model spectra, based on work reported by Lanz et al. (1996) and Barstow, Hubeny 
& Holberg (1998, 1999), were calculated using the non-LTE code TLUSTY (Hubeny & 
Lanz 1995). For this initial analysis, we hxed the stellar temperature and surface gravity 
(Teff = 54, 000 K, logg = 7.5) at the grid points closest to the values determined using the 
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Balmer and Lyman lines (Barstow, Hnbeny & Holberg 1998). Apart from the He abundance, 
which was allowed to vary freely between the grid limits of 10“^ and 10“®, the heavy element 
abundances were hxed at values determined in earlier homogeneous-mo del analyses of G191- 
B2B (C/H = 4.0 X 10-^ N/H = 1. x 10“^ 0/H = 9.6 x 10"^ Si/H = 3.0 x 10“^ P/H = 
2.5 X 10“®, S/H = 3.2 X 10“^ Fe/H = 1.0 x 10-^ Ni/H = 5.0 x 10"^). 

The value taken for Fe/H lies between limits established by FUV (2.4 x 10“®, Venues & 
Lanz (2001)) and EUV (3—4x 10“^, Barstow, Hnbeny & Holberg (1999)) analyses. The effect 
of Fe/H in the 225 — 245 A band is to change the level of the overall spectrum, and we have 
verihed that this does not affect the conclusions reached in our analysis. The ISM H I and He 
I column densities were hxed at values obtained from analysis of the broader-band, lower- 
resolution EUVE spectrum (Barstow, Hnbeny & Holberg 1999); H I = 2.15 x 10^®cm“^, 
He I = 2.18 X lO^^cm”^. The parameters varied during the ht were the column density 
of He H in the line-of-sight (Nneii) and the photospheric He abundance (nne: measured by 
numbers of nuclei). For this ht, the data were summed in bins of width 0.060 A, equivalent 
to R = 4000. Given the uncertainty in the resolving power (section 2, R = 3000 — 4000), the 
hts were performed also for lower values of R, but yielded no signihcant change in the results 
presented below in section 4. The best ht to the data, shown by the red line in Figure 1, was 
obtained for Nneii = 5.97 x lO^^cm”^ and nne = 1-60 x 10“^, and in Figure 2 we show the 
1, 2, 3, 5 and 10 sigma contours for the two parameters. The 3 sigma contour is the locus 
on which exceeds the minimum by 11.8, and within which the parameter conhdence level 
is greater than 99.7 percent. We use this contour to derive 99-percent conhdence limits of 
5.76 - 6.18 X lOi^cm-2 for Nneii and 1.31 - 1.91 x 10"^ for nne- 


4. Discussion 

The good agreement between the best-ht model and the data in Figure 1 is striking, for 
example at the prominent absorption feature at 233.5 A produced by a cluster of O IV lines. 
Many other features are present, which are mainly blends of large numbers of Fe V and Ni 
V lines. The broad features between 227 and 232 A are a characteristic of the overlapping 
series of interstellar He H absorption lines superposed on a continuum. This is shown more 
clearly in the upper panel of Figure 3, an expanded view of the region 226 — 232 A in Figure 
1. Taken with the strong depression of the hux below 227 A this is strong evidence that 
interstellar He H is present along the line of sight. Gonclusive proof is obtained when the 
data are htted by a model in which Nneii is set to zero. The degradation of the ht is evident 
in the lower panel of Figure 3, particularly in the region below 229 A . Further, in this case 
the best-ht value of 8.0 x 10“® for nne is about four times the upper limit we obtained from 
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Wavelength (A) 

Fig. 1.— High resolution EUV spectrum of G191-B2B, obtained with the J-PEX spectrom¬ 
eter in the waveband 221 — 244 A. All data points have error bars. The red histogram is 
the best-£t model of the star and ISM. The strongest predicted lines of He, C, N, O, and 
P are labeled with their ionization state and wavelength. Lines of Fe and Ni, too numer¬ 
ous to include here, account for some unlabelled individual features and broader absorption 
structures. 



Fig. 2.— The 1, 2, 3, 5 and 10 sigma contours for the two parameters varied in fitting 
a uniform abundance white dwarf atmosphere model to the measured spectrum of G191- 
B2B. The 3 sigma contour corresponds to a confidence level of 99.7 percent. The secondary 
contour at 10 sigma locates a weak secondary minimum in the distribution. 
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the absence of detectable He II at 1640 A in the STIS spectrum of G191-B2B. On the other 
hand, the best-£t model in the upper panel of Figure 3, in which UHe = 1-6 10“^, is consistent 
with the STIS limit. 

The best-fit He H column density of 5.97 x lO^^cm”^ implies a He ionization fraction, 
based on EUVE measurements of the He I column density, of ~ 0.73. This is substan¬ 
tially higher than the typical range of 0.25-0.50 in the local interstellar medium (LISM; e.g. 
Barstow et ah 1997). However, a possible circumstellar (CSM) component has been identi¬ 
fied recently through analysis of the CIV absorption line doublet (1548.202 and 1550.774 A) 
in the STIS FUV spectrum of G191-B2B (Bannister et al. 2001). Therefore some of the He 
H detected in the J-PEX spectrum may be circumstellar, although the fraction would have 
to be at least one third to bring the interstellar component within the LISM range. The 
known LISM and CSM components are separated by ~ 8kms“^, and therefore could not be 
resolved by J-PEX. 

The photospheric absorption line at 243.026 A is predicted to be the strongest in the 
J-PEX waveband. Unfortunately this lies in a region of reduced exposure (Section 3), and 
Figure 1 shows no feature in the measured spectrum at this wavelength. However, the depth 
of the predicted line is similar to the magnitude of the statistical errors for the data points, 
yielding only a very weak constraint on the possible He abundance. A weak absorption line 
is seen at the position of He H 237.331 A but likewise the photon statistics do not allow a 
significant detection. 

The above discussion of the column density of ISM/CSM He H, and the evidence for 
photospheric He, is underpinned by the good agreement between the homogeneous stellar 
atmosphere models and the observational data, as shown by the good correspondence be¬ 
tween predicted and observed features in Figure 1 and the upper panel of Figure 3. However, 
closer inspection also reveals several significant features present in the observed but not in 
the synthetic spectrum, for example at 229.3 and 231.2 A . This indicates that the model 
atmosphere is incomplete and that other elements should be included in addition to G, N, 
O, Si, P, S, Fe and Ni. 


5. Conclusions 

We have presented our first analysis of the high-resolution EUV spectrum of G191-B2B 
obtained with the J-PEX spectrometer, which has the highest resolving power (3000-4000) 
achieved so far in the EUV and X-ray wavebands. The results show conclusively that ionised 
He is present along the line-of-sight, and yield a column density of 5.97 x 10^^ cm“^. However 
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Fig. 3.— Upper panel: Expanded view of the J-PEX spectrnm of G191-B2B in the wave¬ 
length range 226 — 232 A spanning the He II Lyman series limit. The red histogram is the 
best-htting model of the stellar atmosphere and absorption in the ISM. 

Lower panel: Expanded view of the J-PEX spectrnm of G191-B2B in the range 226 — 232 
A bnt this time showing the best £t, traced by the green histogram, of a model in which 
absorption by interstellar and circnmstellar He has been removed. 









































































- 9 - 

if this is all interstellar, the measured column density yields an estimated He ionisation 
fraction signihcantly higher than is typical of the LISM. Some of this material may be 
associated with a newly discovered circumstellar component, but further investigation is 
necessary to demonstrate whether this is plausible. 

The white dwarf model which best fitted the data included a significant abundance, 
1.6 X 10“^, of photospheric He, and at the 99-percent confidence level we could exclude 
models containing no photospheric He. However, we caution that a simple model, in which 
elements are distributed uniformly in the atmosphere, was used in this analysis, and further 
work with models in which stratihcation of elements is allowed, and taking greater advantage 
of other G191-B2B observations, is needed to verify or disprove our result. 
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